Limits of quantum speedup in photosynthetic light harvesting 
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It has been suggested that excitation transport in photosynthetic light harvesting complexes features speedups 
analogous to those found in quantum algorithms. Here we compare the dynamics in these light harvesting 
systems to the dynamics of quantum walks, in order to elucidate the limits of such quantum speedups. For the 
Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria, we show that while there is indeed speedup at 
short times, this is short lived (70 fs) despite longer lived (ps) quantum coherence. Remarkably, this time scale 
is independent of the details of the decoherence model. More generally, we show that the distinguishing features 
of light-harvesting complexes not only limit the extent of quantum speedup but also reduce rates of diffusive 
transport. These results suggest that quantum coherent effects in biological systems are optimized for efficiency 
or robustness rather than the more elusive goal of quantum speedup. 



I. INTRODUCTION 

In the initial stages of photosynthesis, energy collected 
from light is transferred across a network of chlorophyll 
molecules to reaction centers 0] [2). Recent experimental 
evidence showing long lived quantum coherences in this en- 
ergy transport in several photosynthetic light-harvesting com- 
plexes suggests that coherence may play an important role 
in the function of these systems fl3}{8l. In particular, it has 
been hypothesized that excitation transport in such systems to 
feature speedups analogous to those found in quantum algo- 
rithms |3][9). These comments have attracted much interest 
from quantum information theorists ll9 ffl"8l . although clearly 
photosynthesis is not implementing unitary quantum search 
191 . The most direct analogy to such transport is found in 
quantum walks, which form the basis of a powerful class of 
quantum algorithms including quantum search [19-24]. Un- 
like idealized quantum walks, however, real light harvesting 
complexes are characterized by disorder, energy funnels and 
decoherence. Whether any quantum speedup can be found in 
this situation has remained unclear. 

Quantum walks are an important tool for quantum algo- 
rithms [19-24]. On the line, quantum walks feature ballis- 
tic spreading, (x 2 ) oc t 2 , compared to the diffusive spreading 
of a classical random walk, (x 2 ) oc t, where (x 2 ) denotes 
the mean squared displacement. Moreover, because they use 
superpositions instead of classical mixtures of states, on any 
graph with enough symmetry to be mapped to a line, quantum 
walks spread along that line in linear time - even when classi- 
cal spreading is exponentially slow [20]. We shall refer to this 
enhanced rate of spreading as a generic indicator of "quantum 
speedup." Such quantum speedup is important in quantum 
information processing, where it can lead to improved scal- 
ing of quantum algorithms relative to their classical alterna- 
tives. Examples of such algorithmic speedup include spatially 
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structured search ETt element distinctness ll22l and evaluat- 
ing AND-OR formulas ll23l . Quantum walks even provide a 
universal implementation for quantum computing l24l . 

Quantum walks also constitute one of the simplest models 
for quantum transport on arbitrary graphs. As such they pro- 
vide a theoretical framework for several physical processes, 
including the transfer of electronic excitations in photosyn- 
thetic light harvesting complexes [3j[9)- The closed system 
dynamics in a light-harvesting complex are generally well de- 
scribed by a tight-binding Hamiltonian that is restricted to the 
single excitation subspace [2 |, 

H = J2 E n\n)(n\ + ^ Jnm\n)(m\. (1) 

Here \n) represents the state where the nth chromophore (site) 
is in its electronic excited state and all other chromophores 
are in their ground state. E n is the electronic transition en- 
ergy of chromophore n and J mn is the dipole-dipole coupling 
between chromophores n and m. This is a more general vari- 
ant of the Hamiltonian for standard continuous-time quantum 
walks, where E n = and J mn g {0, 1}. The salient differ- 
ences of the general model, variable site energies and cou- 
plings, arise naturally from the structure and role of light- 
harvesting complexes. Non-constant site energies can serve 
as energy funnels and can enable the complex to absorb at 
a broader range of frequencies, while variable couplings be- 
tween sites reflect their physical origin as dipole-dipole inter- 
actions. These differences yield excitation dynamics that can 
deviate significantly from quantum walks. 

Here we study the key question of whether excitation trans- 
port on light-harvesting complexes shows quantum speedup. 
Such quantum speedup would be necessary for any quan- 
tum algorithm that offers algorithmic speedup relative to a 
classical search of physical space. (Achieving true algorith- 
mic speedup, i.e., improvement over the best classical al- 
gorithms, would also require suitable scaling of the space 
requirements 1251 .") We address this with a study of the 
Fenna-Matthews-Olson (FMO) complex of green sulfur bac- 
teria, a small and very well-characterized photosynthetic com- 
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FIG. 1: (a) Crystal structure of FMO complex of C. tepidum (Protein 
Data Bank accession 3ENI), with lines between the chromophores 
representing dipolar couplings. The thickness of the lines indicates 
the coupling strengths. Only couplings above 15 cm -1 are shown; 
the largest coupling is 96 cm -1 . The full Hamiltonian is given in[X| 
(b) Site energies, shown relative to 12 210 cm - , in the reduced di- 
mensionality model derived from mapping the strongest couplings 
onto a one-dimensional graph. The red sites (1,6) are source sites 
at which the excitation enters the complex and the black site (3) is 
the trap site from which the excitation is transferred to the reaction 
center 12711281 . 

plex JT] 12 |26) l27l . the same complex for which long lived 
quantum coherences were recently observed and suggested to 
reflect execution of a natural quantum search Q. We also 
consider the dynamics of transport along a theoretical model 
of an extended chain of chromophores to elucidate the sys- 
tematic influence of variable site energies and dephasing on 
transport. 



H. FENNA-MATTHEWS-OLSON COMPLEX 

The FMO complex acts as a quantum wire, transporting 
excitations from a large, disordered antennae complex to a 
reaction center. In addition to the possibility of quantum 
speedup (3j |9), recent studies have speculated that coher- 
ence may assist unidirectional transport along this wire [5| 
or suggested that it may contribute to overall efficiency ATI . 
The crystal structure of FMO shows three identical subunits 
that are believed to function independently, each with seven 
bacteriochlorophyll-a molecules embedded in a dynamic pro- 
tein cage l26l . A refined model Hamiltonian for the single 
excitation subspace is available from detailed quantum chem- 
ical calculations |27|, and the orientation of the complex was 
recently verified experimentally [28 1. By neglecting the weak- 
est couplings in this model Hamiltonian, we see that transport 
in an individual monomer of FMO can be mapped to a one- 
dimensional path between chromophores, as shown in Fig- 
ure [TJ This mapping of the excitation transport in FMO to a 
single dimension allows us to make contact with known re- 
sults for quantum transport in one dimension and to quantita- 
tively assess the extent of quantum speedup. 



III. COHERENT DYNAMICS 

Under Hamiltonian dynamics, quantum walks on highly 
symmetric graphs can spread ballistically, but any signifi- 



cant lack of symmetry can lead to localization. Consider 
transport along the infinite line. Here, a random variation 
of any magnitude in the site energies leads to Anderson lo- 
calization (29]. Similarly, random variations in the coupling 
strengths contribute to localization l30l . Systematic variations 
in site energies or coupling strengths can also cause localiza- 
tion ||3T| , with the well-known instance of Bloch oscillations 
and Stark localization deriving from a linear bias in site ener- 
gies P2| . Both systematic and random variations in site en- 
ergies remove the symmetry necessary for Bloch's theorem, 
so it is not surprising that their combination leads to localiza- 
tion as well l33l l34l . Adding disorder into the graph struc- 
ture also usually causes localization [ 35 1 (although such dis- 
order can also reduce localization if it makes an already dis- 
ordered graph more connected [36]). The varied couplings, 
energy funnel and disordered energies evident in the FMO 
Hamiltonian depicted in Figure [TJ suggest that localization 
due to several of these effects will be significant for light- 
harvesting complexes. The standard measure of localization 
under coherent dynamics, the inverse participation ratio l29l 
£ = Y^,i IV^I 2 / Si IV'il 4 ' f° r tne amplitudes ipi in the site ba- 
sis of an eigenstate ip of our model Hamiltonian confirms this 
intuition, since the typical eigenstate for FMO occupies only 
£ - 2 sites (see also Ref. El). 

We can estimate the timescale for localization t\ oc based 
upon the experimentally accessible parameters of geometry, 
localization length £ and an average coupling strength J. 
Consider the speed of the quantum walk as an upper bound 
on excitation transport speed prior to localization. For the 
continuous -time quantum walk on the infinite line, it is well 
known that (x 2 ) ~ 2J 2 t 2 /h 2 |38|. Similarly, starting at one 
end of an infinite line, (x 2 ) ~ 3J 2 t 2 /h 2 . This gives an av- 
erage speed gJ/h, where g = y2 and g = y3 respectively, 
yielding the bound ti oc > h£/gj. Our simulations of strict 
Hamiltonian dynamics (no decoherence) with FMO as de- 
scribed below find the onset of localization at ii oc ~ 70 fs, 
which is close to the bound h£/gj ~ 100 fs from the mean in- 
verse participation ratio £ avg w 2, the mean coupling strength 
J avg ~ 60cm -1 and g — y3 (this value for g is suggested 
by the dominant pathways in Figure[JJfor transport starting at 
sites 1 or 6, which are the chromophores closest to the antenna 
complex 12711281 "). 



IV. DECOHERENCE 

Decoherence, i.e. non-unitary quantum evolution, is also an 
essential feature of excitation dynamics in real systems. In 
light harvesting complexes decoherence arises from interac- 
tions with the protein cage, the reaction center and the sur- 
rounding environment. Recently it has been shown that some 
degree of fluctuation of electronic transition energies (i.e. de- 
phasing in the site basis) increases the efficiency of transport 
in FMO and other simple models otherwise limited by Ander- 
son localization. The intuition is that at low levels dephasing 
allows escape from localization by removing destructive in- 
terference, but at high levels it inhibits transport by inducing 
the quantum Zeno effect l9l [lOl [Til [131 . 
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We consider here two essential types of decoherence: de- 
phasing in the site basis and loss of excitations. Under the 
Born-Markov approximation, time evolution with this deco- 
herence model is determined by the system Hamiltonian H 
and a sum of Lindblad operators C according to 1391 



(2) 



^ = -l[H ) p]+Ci 0SS (p) + C deph (p). 



Class describes loss of excitations with site dependent rates 7„, 
including trapping to a reaction center, and is specified by the 
operator 

(n|£i 0S8 (p)|m) = - 7 " - 7m (n\p\m), (3) 

which results in exponential decay of the diagonal elements of 
the density matrix p (and corresponding decay of off-diagonal 
elements, which ensures positivity). £deph describes dephas- 
ing, which is specified by the operator 



-(1-S nm )(n\p\m) (4) 



(n|£deph(p)M 



with site dependent dephasing rates T„, and yields exponen- 
tial decay of the off -diagonal elements of p. 



A remarkable feature of our results for FMO is that, as we 
shall demonstrate, they are independent of the finer details of 
the bath dynamics and system-bath coupling. Thus for sim- 
plicity, we specialize here to the Haken-Strobl model |40| 
and restrict losses to trapping by the reaction center. The 
ease of calculations with the Haken-Strobl model has made it 
popular for simulating the dynamics of light harvesting com- 
plexes (lOlEliniSIlHa. Following Rebentrost et al. ifTol . 
we use the spatially uniform, temperature dependent dephas- 
ing rate T — 2irkBTEn/h 2 uj c , for an ohmic spectral den- 
sity with bath reorganization energy En = 35 cm -1 and cut- 
off frequency uj c = 150 cm" 1 . This dephasing rate holds 
when chromophores are treated as qubits coupled to indepen- 
dent reservoirs in the spin-boson model, in the thermal (or 
Markovian) regime t h/ksT 11391 . which corresponds to 
t ^> 25 fs at 300 K. Although the exact parameters of the 
protein environment surrounding FMO are unclear (5), these 
are reasonable estimates. This model gives dephasing rates 
(69 fs)" 1 at 77 K and (18 fs)- 1 at 300 K. Also, as in Ref. ED 
we restrict trapping to site 3 at a rate of 73 = 1 ps _1 . Finally, 
we neglect exciton recombination since it occurs on much 
slower timescales (^1 ns -1 ) and in any case it does not alter 
the mean-squared displacement. Therefore, we set 7^ = for 
i 7^ 3. While more realistic decoherence models would incor- 
porate thermal relaxation, spatial and temporal correlations in 
the bath and strong system-bath coupling [43 1, this treatment 
is sufficient for analyzing the quantum speedup, as justified 
in detail by comparison with the corresponding analysis using 
more realistic simulations of FMO [43 1 in[B] 



V. LIMITS OF QUANTUM SPEEDUP 

To access the extent of quantum speedup for a system 
that can be mapped to a line, a natural measure is the expo- 
nent b of the power law for the mean squared displacement 
(x 2 ) oc t b . In particular, we use the best-fit exponent b from 
the slope of the log-log plot of the mean squared displace- 
ment (x 2 ) = Ti'[px 2 ]/ Tr p versus t. The value b = 1 cor- 
responds to the limit of diffusive transport, whereas b = 2 
corresponds to ideal quantum speedup as in a quantum walk 
(ballistic transport). 

The timescale for quantum speedup is generally bounded 
above by both the timescale for dephasing, ideph = 1 /T, and 
the timescale for static disorder to cause localization, t] oc . To 
see this, it is illustrative to consider transport along a linear 
chain with constant nearest neighbor couplings J. For dephas- 
ing but no static disorder or energy gradient (E„ = 0), the 
mean squared displacement shows a smooth transition from 
ballistic to diffusive transport BUBS 1 > 



4J 2 

n 2 r 



1 , 



,-rt 



(5) 



The corresponding power law transition is shown in Figure|2^. 
Technically, transport remains super-diffusive even after ideph, 
since the power never drops below 6=1. In contrast, for 
static disorder but no dephasing, there is a sudden transition 
from ballistic transport to essentially no transport at all (local- 
ization). The power law starts at b = 2, then drops and begins 
to oscillate wildly after t\ oc as the wave function continues to 
evolve in a confined region. Several examples are shown in 
Figure for variable strengths of disorder. Figure [2J; shows 
the behavior with both strong static disorder and dephasing, 
as in light harvesting complexes. In this case, the transport 
can even exhibit a sub-diffusive power law. The reason for 
this sub-diffusive behavior may be easiest to understand by 
analogy to the Anderson model (random site energies) in an 
infinite chain. Consider transport in such a system under weak 
dephasing and with an ensemble average over different real- 
izations of strong static disorder. In this case, one expects 
transport must transition from ballistic (b = 2) at short times, 
to localized (b = 0) at intermediate times, to diffusive at long 
times (b = 1). (The situation at long times is analyzed more 
explicitly in Section VI and [C]) For stronger dephasing as in 
Figure |2j; (and as we shall see with FMO), the fully localized 
regime with no transport may never be realized, but transport 
will still be sub-diffusive for intermediate times. 

Figure [3] shows the results of a simulation of FMO dy- 
namics made with our simple decoherence model. The dis- 
placement Xi of a site i is given by its position in the one- 
dimensional mapping that is presented in Figure [TJ). Note that 
our simulations use the full Hamiltonian; the one-dimensional 



1 This result, obtained in the context of excitation transport in molecular 
crystals, also applies to decoherent continuous-time quantum walks since 
they use the same tight-binding Hamiltonian. 
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FIG. 2: The best fit power law exponent b for the mean-squared dis- 
placement with J/h = 1, for transport along a linear chain with 
either (a) increasing dephasing (T = 1, 3, 9) and no static disorder 
(E„ = 0), (b) no dephasing (r = 0) and increasing static disorder 
(independent site energies E n from a single Gaussian distribution 
normalized to standard deviation a = 1, 3, 9) or (c) same disorder as 
(b), but with finite dephasing rate T = 1. Panels (b) and (c) show 
results for single instances of disorder. When an ensemble average 
over different realizations of static disorder is taken, the power law 
varies smoothly instead of oscillating as in (b) and (c). 



map is only used in the analysis of the results, to determine the 
displacement xi of each site. The upper-left panel, (a), shows 
the mean squared displacement (a; 2 ) of an excitation initially 
localized on site 6, one of the two excitation source sites, and 
the lower-left panel, (b), the best fit exponent b for the power 



law 



t . Our results show that even though coher- 



ence in this simple model lasts for ~ 500 fs (see Figure |3};), 
a transition from initially ballistic to sub-diffusive transport 
occurs after only ~ 70 fs, independent of the dephasing rate. 
While increased dephasing causes a faster initial decrease of 
the power law, dephasing alone cannot lead to sub-diffusive 
transport, as is clear from Equation Q. We note that the tran- 
sition to sub-diffusive transport at 70 fs occurs at the same 
time-scale as the onset of localization (Section III i, imply- 



ing that even though there is persistent coherence beyond this 
time, it no longer yields quantum speedup because of static 



disorder. The sub-diffusive power law at intermediate times 
(~ 100 fs - 2ps) arises from the interplay of this disorder- 
induced localization and dephasing that is discussed above for 
the linear chain model. At longer times the power law expo- 
nent b goes to zero because of the finite size of the system. 
Since complete energy transfer through FMO takes picosec- 
onds, this analysis shows that most of the excitation transport 
is formally sub-diffusive. 

While FMO is a relatively small system, so that terms such 
as "ballistic" and "diffusive" cannot literally describe trans- 
port across its seven chromophores, this time dependence of 
the power law of spreading would also characterize larger ar- 
tificial or natural systems. Our results here are robust to vari- 
ations in the strength of the trap, and to whether the initial 
excitation is at site 1 or 6, the sites believed to be the primary 
source for excitations in FMO {271 . They are also indepen- 
dent to variations of the map to a one-dimensional system in 
Figure [T] such as counting the number of chromophores away 
from the trap site instead of along the entire line, or using the 
real space distance between adjacent chromophores instead of 
assuming that each pair is merely separated by a unit lattice 
distance. 

Our timescale for short-lived quantum speedup in FMO 
should be largely independent of the details of the decoher- 
ence model, since the non-uniform energy landscape will al- 
ways limit ballistic transport to times prior to ti oc ~ 70 fs. 
This is confirmed by applying our analysis to the results of 
recent calculations for FMO made with a considerably more 
sophisticated decoherence model that incorporates thermal re- 
laxation, temporal correlations in the bath and strong system- 
bath coupling l43l . As demonstrated in|B] these more realistic 
calculations also predict a loss of quantum speedup after 70 fs. 



VI. DIFFUSIVE TRANSPORT 

To gain general insight into the interplay between disorder 
and dephasing for light harvesting complexes, we now con- 
sider the dynamics of an infinite linear chain at long times 
with variable site energies, couplings and dephasing rates. 
By extending an analysis for damped Bloch oscillations (46), 
we find that any set of non-zero dephasing rates T n asymp- 
totically leads to diffusive transport (x 2 ) ~ 2Dt. Haken 
and Reineker also performed a reduction of the Haken-Strobl 
model without disorder to a diffusion equation along similar 
lines l47l . Note that since the Haken-Strobl model does not 
include energy relaxation, no classical drift velocity will be 
obtained from these dynamics. The time evolution of an arbi- 
trary density matrix element for an infinite linear chain under 
site dependent dephasing rates is given by 



±\Jn—\Pn—\,m JnPn+l,m 

JmPn,m+l Jm— lPn,m— l) 

i r„ + r 

■^(E n ~E m )-\ (1 — 5 nm ) p n ,m, (6) 
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FIG. 3: Results of simulations on FMO with the decoherence model described in the text, (a) Log-log plot of mean-squared displacement (x 2 ) 
as a function of time, for an initial excitation at site 6. (b) Power law exponent for the mean-squared displacement, given by the slope of the 
plot in panel a. The dashed line 6=1 separates super- and sub-diffusive transport (see text), (c) Total site coherence, given by the sum of the 
absolute value of all off-diagonal elements of the density matrix in the site basis, indicates persistent coherence for ~ 500 fs. (d) Oscillating 
site populations (shown for 77 K only) also indicate persistent coherence for hundreds of femtoseconds. 



where J„ = J n ,„+i- If we neglect terms in the second off- 
diagonal relative to the diagonal terms (see justification in[C]> 
we obtain 

Pn+l.n — —^Jn(Pn,n ~ Pn+l,n+l) ~ ^n)Pn+l,n, 

(7) 

where A„ = E n+1 - E n and 1^ = (T n + T n+1 )/2. This 
equation has the solution 



Pn+l,n — T e 

n 



-(iA„/ft+r n )t 



x / J„(p„ + i,„+i - p n ,n)e {iAn/H+rn)t 'dt' (8) 



,Jn\Pn+l,n+l ~ Pn,n) 

'' A„ - ihY' n ' 



(9) 



in the nearly stationary regime where site populations vary 
slowly compared to 1/r^ and H/A n . Inserting (|9]l and the 
corresponding result for p n -x,n m lo ® we obtain 



This is now a classical random walk with variable bond 
strengths between sites. For asymptotically large times it has 
been proven to approximate the diffusion equation with coef- 
ficient given by 



D = 




(11) 



as long as the average over sites in the diffusion coefficient is 
well-defined |48|. With A„, J„ and r„ constant, this matches 
the known diffusion coefficient ll46l l49l . For constant cou- 
pling J n = J and constant dephasing T' n = T, we note that 
( fTTj ) reduces to 



D = 



2J 2 T 



(12) 



which is plotted as a function of (A^ ) 1 / 2 and T in Figure [4] 
2 . Note that with (A 2 ) = 0, this is the long time limit of 
(|5). In the case of Equation ( fl2| ), for a given degree of static 



2 7 2 r' 

Pnn — , 2 fc2p/2 \Pn+l,n+l Pnn) 



9 F' 

+ 1 fi2r/2 (Pn-l.n-1 ~ Pnn)- (10) 

^n-l ' n 1 n-1 



2 We note that the parameter (A 2 ) 1 / 2 does not completely specify the de- 
gree of localization, although it does show that Stark and Anderson local- 
ization equivalently influence diffusive transport. For example, the inverse 
participation ratio in these two cases differs, even with the same value of 
<A 2 >V 2 . 
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FIG. 4: Diffusion coefficients at long times as a function of a constant 
dephasing rate V and the static disorder between adjacent sites, as 
calculated by {12) with J = h = 1. Solid lines are contours; the 
dashed line is the optimal dephasing rate for given static disorder. 
The ideal quantum walk, which is non-diffusive, is at the origin. 

disorder, we find an optimal dephasing rate hT — (A 2 ) 1 / 2 . 
However, Figure [4] shows that for a given dephasing rate, ex- 
citation transport would be improved by reducing the static 
disorder to further delocalize the system. 



VII. CONCLUSIONS 

These principles show that quantum speedup across a pho- 
tosynthetic or engineered system requires not only long-lived 
quantum coherences but also excitons delocalized over the en- 
tire complex. Such completely delocalized excitons do not 
exist in FMO and are also unlikely in other light harvesting 
complexes because of the presence of energy gradients and 
disorder. Moreover, Equation (JT2j makes it clear that the con- 
ditions needed for longer lived quantum speedup (reduced de- 
phasing and static disorder) are those necessary for faster dif- 
fusive transport as well. 

The short-lived nature of quantum speedup in light harvest- 
ing complexes that we have established here implies that the 
natural process of energy transfer across these complexes does 
not correspond to a quantum search. Indeed, neither a formal 
quantum nor classical search may be necessary, since pre- 
determined (evolved) energy gradients can guide relaxation 
to reaction centers, even though such gradients suppress co- 
herent and dephasing-assisted transport. This is particularly 
relevant for systems such as FMO which either receive exci- 
tations one at a time or have isolated reaction centers. Instead 
of yielding dynamical speedup like that in quantum walk algo- 
rithms, quantum coherence in photosynthetic light harvesting 
appears more likely to contribute to other aspects of transport, 
such as overall efficiency or robustness. We emphasize that 
a restricted extent of quantum speedup does not imply that 
there is no significant quantum advantage due to long-lived 
coherence in electronic excitation energy transfer. Identifying 
the specific nature of any "quantum advantage" for FMO will 
clearly require more detailed analysis of the dynamics, partic- 
ularly in the sub-diffusive regime. Related examples of such 



"quantum advantage" are found in the LH2 complex of pur- 
ple bacteria, where coherence has been specifically shown to 
improve both the speed ll50l and robustness |51 1 of transport 
from the B800 to B850 ring. 
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Note added 

While revising this manuscript, a reduction of the Haken- 
Strobl model to a classical random walk similar to the one we 
perform in Section[Vl]was published by Cao and Silbey ll52l . 



Appendix A: FMO Hamiltonian 



In all of our calculations, we use the Hamiltonian calculated 
for C. tepidum by Adolphs and Renger 11271 . In the site basis, 
this is given by 
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with units of cm 1 and a total offset of 12 210 cm 1 . Bold 
entries indicate those shown in Figure [T] In units with H = 1, 
we note that the rate 1 ps _1 = 5.3 cm . 



Appendix B: Reduced hierarchy equations model for FMO 

Energy transfer dynamics in photosynthetic complexes can 
be difficult to model because perturbations from the surround- 
ing protein environment can be large, and the timescale of 
the protein dynamics is similar to the timescales of excitation 
transport. This makes common approximations involving per- 
turbative treatments of system-bath coupling and Markovian 
assumptions on the bath invalid. Recently a non-perturbative, 
non-Markovian treatment of energy transfer has been formu- 
lated by Ishizaki and Fleming |43 1. This model assumes: (i) a 
bilinear exciton-phonon coupling, (ii) protein fluctuations that 
are described by Gaussian processes, (iii) a factorizable initial 
state of chromophores and protein environment, (iv) protein 
fluctuations that are exponentially correlated in time, and (v) 
no spatial correlations of fluctuations. 
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FIG. 5: For the reduced hierarchy equation model as applied to FMO, 
we plot the power law for the mean squared displacement as in Fig- 
ure [3J). The plot is interpolated from a fourth order polynomial fit of 
density matrix elements for results from a numerical simulation on 
FMO with sampling approximately every 8 fs, with the initial excita- 
tion at site 6 0. 

To verify our results in Figure [3] of the main paper with 
this more realistic model for excitation dynamics, we apply 
the power law analysis of mean squared displacement to the 
results of the simulations that were used to calculate entangle- 
ment dynamics in Ref. |fl6l using the non-perturbative non- 
Markovian method of Ref. |5). These simulations used a re- 
organization energy of the protein environment of 35 cm -1 , 
phonon relaxation time of 100 fs, and a reaction center trap- 
ping rate of (4 ps)" 1 , all of which are consistent with litera- 
ture on FMO Il27ll53i . Figure|5]shows the power law analysis 
for 2 different temperatures, 77 K and 300 K. Our analysis of 
these realistic simulations show that longer lasting coherences 
(as discussed in Refs. 15]|43]) are evident in the power law os- 
cillations, but that transport is still nevertheless sub-diffusive 
after ~ 70 fs, for both temperatures. 



It is particularly noteworthy that our results for FMO hold 
even under a model incorporating finite temperature relax- 
ation. In principle, thermal relaxation could bias the dynamics 
on a linear chain with a classical drift velocity |54| to give 
the power law (a; 2 ) oc t 2 even without quantum speedup. 
However, the above analysis of FMO simulations made with 
the most realistic treatment of relaxation available today (Fig- 
ure [5]) shows that this is not the case in this system. 



Appendix C: Decay of coherences in linear transport 

For a constant dephasing rate T ^> J/h, there is a sim- 
ple analytical argument that the off-diagonal elements of the 
density matrix in the site basis can be neglected relative to the 
main diagonal l55ll . In the case of transport along a line, appli- 
cation of the analysis of Ref. Il55ll to Equation (|6]i shows that 
the first off-diagonal elements decay at rate T, unless more 
coherence is generated from the main diagonal. If the main 
diagonal elements have order 1, then the first off-diagonal el- 
ements cannot have magnitude greater than order J/hT. A 
similar argument bounds the second diagonal elements as less 
than {J/HT) 2 and so forth. For large T, this justifies neglect- 
ing this second off-diagonal relative to the main diagonal. The 
higher order terms in this expansion can also be calculated ex- 
plicitly 11521 . 

From extensive numerical experiments with infinite chains, 
we have found uniformly that the off-diagonal elements of 
the density matrix decay nearly monotonically after time 1 /T 
even when Y ^ J/h does not hold, although an analytic proof 
of this result has eluded us. Several examples are presented in 
Figure [6] This matches the analytical result for finite systems 
in the Haken-Strobl model that coherence must vanish even- 
tually for any non-zero Y |56|. 



[1] R. E. Blankenship, Molecular mechanisms of photosynthesis 

(Wiley-Blackwell, Maiden, MA, USA, 2002). 
[2] Y.-C. Cheng and G. R. Fleming, Dynamics of Light Harvesting 

in Photosynthesis, Annu. Rev. Phys. Chem. 60, 241 (2009). 
[3] G. S. Engel et ai, Evidence for wavelike energy transfer 

through quantum coherence in photosynthetic systems, Nature 

446, 782 (2007). 

[4] H. Lee, Y.-C. Cheng, and G. R. Fleming, Coherence Dynamics 
in Photosynthesis: Protein Protection of Excitonic Coherence, 
Science 316, 1462 (2007). 

[5] A. Ishizaki and G. R. Fleming, Theoretical examination of 
quantum coherence in a photosynthetic system at physiological 
temperature, Proc. Natl. Acad. Sci. USA 106, 17255 (2009). 

[6] T. R. Calhoun et al, Quantum Coherent Enabled Determina- 
tion of the Energy Landscape in Light-Harvesting Complex II, 
J. Phys. Chem. B. Lett. 113, 16291 (2009). 

[7] E. Collini et al. , Coherently wired light-harvesting in photosyn- 
thetic marine algae at ambient temperature, Nature 463, 644 
(2010). 

[8] G. Panitchayangkoon et al., Long-lived quantum coherence in 



photosynthetic complexes at physiological temperature, (2010), 
larXiv:1001 5108vl| 
[9] M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru- 
Guzik, Environment-Assisted Quantum Walks in Photosyn- 
thetic Energy Tra nsfer, J. Chem. Phys. 129, 174106 (2008), 
larXiv:0805.274"T1 

[10] P. Rebentrost, M. Mohseni, I. Kassal, S. Lloyd, and A. Aspuru- 
Guzik, Environment-Assisted Quantum Transport, New J. 
Phys. 11, 033003 (2009), |arXlv:0807.0929| 

[11] P. Rebentrost, M. Mohseni, and A. Aspuru-Guzik, Role of 
quantum coherence in chromophoric energy transport, J. Phys. 
Chem. B 113 (2009), |a7^GvT08 06.4725 

[12] M. Plenio and S. Huelga, Dephasing assisted transport: Quan- 
tum networks and biomolecules, New J. Phys. 10, 113019 
(2008), [mXiv:0807.4902| 

[13] F. Caruso, A. Chin, A. Datta, S. Huelga, and M. Plenio, Highly 
efficient energy excitation transfer in light-harvesting com- 
plexes: The fundamental role of noise-assisted transport, J. 
Chem. Phys. 131, 105106 (2009), |arXiv:090 1 ,4454| 

[14] F. Caruso, A. Chin, A. Datta, S. Huelga, and M. Plenio, En- 



8 




0.01 0.1 1 10 0.01 0.1 1 10 0.01 0.1 1 10 

Time Time Time 



FIG. 6: Decay of coherences. Log-log plot of the sum of the absolute values fcth off-diagonal elements of the density matrix over time for 
constant dephasing and coupling strengths T — J/h = 1. Panels show the quantum walk with dephasing (a), and Stark (b) and Anderson (c) 
localization with dephasing and (A 2 l } 1//2 = n/2. Note that the main diagonal always sums to 1. 



tanglement and entangling power of the dynamics in light- 
harvesting complexes, (2009), arXiv:0912.0 122| 

[15] A. W. Chin, A. Datta, F. Caruso, S. F. Huelga, and M. B. Plenio, 
Noise-assisted energy transfer in quantum networks and light- 
harvesting complexes, (2009), arXiv:0910.4153 

[16] M. Sarovar, A. Ishizaki, G. R. Fleming, and K. B. Whaley, 
Quantum entanglement in photosynthetic light harvesting com- 
plexes, (2009), |arXiv:0905.3787| 

[17] M. M. Wilde, J. M. McCracken, and A. Mizel, Could light har- 
vesting complexes exhibit non-classical effects at room temper- 
ature?, Proc. R. Soc. A: Math. Phys. (2009), arXiv: 0911.1097| 

[18] K. Bradler, M. M. Wilde, S. Vinjanampathy, and D. B. Uskov, 
Identifying the quantum correlations in light-harvesting com- 
plexes, (2009), arXiv:0912 5ll2vT) 

[19] A. Ambainis, Quantum walks and their algorithmic appli- 
cations, Int. J. Quantum Inf. 1, 507 (2003), arXiv:quant- 

| ph/040 3120 

[20] A. M. Childs et al., Exponential algorithmic speedup by quan- 
tum walk, in 35th Annual ACM Symposium on Theory of 
Computing, pp. 59-68, New York, NY, USA, 2003, ACM, 
|arXiv:qu ant-ph/020913T1 

[21] N. Shenvi, J. Kempe, and K. B. Whaley, Quantum random- 
walk search a lgorithm, Phys. Rev. A 67, 052307 (2003), 
|arXiv:quant-ph/0210 064 

[22] A. Ambainis, Quantum Walk Algorithm for Element Distinct- 
ness, SIAM J. Comput. 37, 210 (2007). 

[23] A. Ambainis, A. M. Childs, B. W. Reichardt, R. Spalek, and 
S. Zhang, Any AND-OR Formula of Size N can be Evaluated in 
time TV 1 / 24 on a Quantum Computer, in 48th Annual IEEE 
Symposium on Foundation of Computer Science, pp. 363-372, 
Los Alamitos, CA, 2007, IEEE Computer Society. 

[24] A. M. Childs, Universal computation by quantum walk, Phys. 
Rev. Lett. 102, 180501 (2009), arXiv:08067l972] 

[25] R. Blume-Kohout, C. M. Caves, and I. H. Deutsch, Climb- 
ing Mount Scalable: Physical Resource Requirements for Scal- 
able Quantum Computer, Found. Phys. 32, 1641 (2002), 
|arXiv:qua nt-ph/0204157l 

[26] R. E. Fenna and B. W. Matthews, Chlorophyll arrangement in a 
bacteriochlorophyll protein from Chlorobium limicola, Nature 
258, 573 (1975). 

[27] J. Adolphs and T. Renger, How proteins trigger excitation en- 
ergy transfer in the FMO complex of green sulfur bacteria, Bio- 



phys. J. 91, 2778 (2006). 

[28] J. Wen, H. Zhang, M. L. Gross, and R. E. Blankenship, Mem- 
brane orientation of the FMO antenna protein from Chlorobac- 
ulum tepidum as determined by mass spectrometry-based foot- 
printing, Proc. Natl. Acad. Sci. USA 106, 6134 (2009). 

[29] P. Phillips, Anderson localization and its exceptions, Annu. 
Rev. Phys. Chem. 44, 115 (1993). 

[30] H. Fidder, J. Knoester, and D. A. Wiersma, Optical proper- 
ties of disordered molecular aggregates: A numerical study, J. 
Chem. Phys. 95, 7880 (1991). 

[31] S. I. Serdyukova and B. N. Zakhariev, Peculiarities of the quan- 
tum interchannel motion: Why the quantum lattice wave does 
not slide down the potential slope, Phys. Rev. A 46, 58 (1992). 

[32] T. Hartmann, F. Keck, H. J. Korsch, and S. Mossmann, Dynam- 
ics ofBloch oscillations, New J. Phys. 6, 2 (2004). 

[33] M. Luban and J. H. Luscombe, Localization by electric fields in 
one-dimensional tight-binding systems, Phys. Rev. B 34, 3674 
(1986). 

[34] A. R. Kolovsky, Interplay between Anderson and Stark Local- 
ization in 2D Lattices, Phys. Rev. Lett. 101, 190602 (2008). 

[35] O. Miilken, V. Pernice, and A. Blumen, Quantum transport 
on small-world networks: A continuous -time quantum walk ap- 
proach, Phys. Rev. E 76, 051125 (2007), |arXiv:0705.1608| 

[36] G. Giraud, B. Georgeot, and D. L. Shepelyanksy, Quantum 
computing of derealization in small-world networks, Phys. 
Rev. E 72, 036203 (2005), |arXiv:quant-ph/0503188l 

[37] T. Brixner et al., Two-dimensional spectroscopy of electronic 
couplings in photosynthesis, Nature 343, 625 (2005). 

[38] N. Konno, Limit theorem for continuous-time quantum walk 
on the line, Phys. Rev. E 72, 026113 (2005), |arXlv7 quant- 
ph/0408T40| 

[39] H.-P. Breuer and F. Petruccione, The theory of open quantum 
systems (Oxford University Press, Oxford, UK, 2002). 

[40] H. Haken and G. Strobl, Exact treatment of coherent and in- 
coherent triplet exciton migration, in The triplet state: Pro- 
ceedings of an international symposium held at the American 
University of Beirut, Lebanon, edited by A. B. Zahlan, pp. 311- 
315, Cambridge University Press, 1967. 

[41] J. A. Leegwater, Coherent versus Incoherent Energy Transfer 
and Trapping in Photosynthetic Antenna Complexes, J. Phys. 
Chem. 100, 14403 (1996). 

[42] K. M. Gaab and C. J. Bardeen, The effects of connectivity, 



9 



coherence, and trapping on energy transfer in simple light- 
harvesting systems studied using the Haken-Strobl model with 
diagonal disorder, J. Chem. Phys. 121, 7813 (2004). 

[43] A. Ishizaki and G. R. Fleming, Unified treatment of quantum 
coherent and incoherent hopping dynamics in electronic energy 
transfer: Reduced heirarchy equation approach, J. Chem. Phys. 
130,234111 (2009). 

[44] E. Schwarzer and H. Haken, The moments of coupled coherent 
and incoherent motion of excitons, Phys. Lett. 42A, 317 (1972). 

[45] M. Grover and R. Silbey, Exciton Migration in Molecular Crys- 
tals, J. Chem. Phys. 34, 4843 (1971). 

[46] A. Kolovsky and A. Ponomarev, Damped Bloch oscillations of 
cold atoms in optical lattices, Phys. Rev. A 66, 053405 (2002), 
|arXiv:qua nt-ph/0206 1 08 

[47] H. Haken and P. Reineker, The Coupled Coherent and Incoher- 
ent Motion of Excitons and its Influence on the Line Shape of 
an Optical Absorption, Z. Physik249, 253 (1972). 

[48] B. D. Hughes, Random Walks and Random Environments Vol. 
2: Random Environments (Oxford University Press, Oxford, 
UK, 1996). 

[49] D. H. Dunlap and V. M. Kenkre, Effect of scattering of the 



dynamic localization of a particle in a time-dependent electric 
field, Phys. Rev. B 37, 6622 (1988). 

[50] S. Jang, M. D. Newton, and R. J. Silbey, Multichromophoric 
Forster Resonance Energy Transfer, Phys. Rev. Lett. 92 (2004). 

[51] Y. C. Cheng and R. J. Silbey, Coherence in the B800 Ring of 
Purple Bacteria LH2, Phys. Rev. Lett. 96 (2006). 

[52] J. Cao and R. J. Silbey, Optimization of Exciton Trapping in En- 
ergy Transfer Processes, J. Phys. Chem. A 113, 13825 (2009). 

[53] M. Cho, H. M. Vaswani, T. Brixner, J. Stegner, and G. R. Flem- 
ing, Exciton analysis in 2D electronic spectroscopy, J. Phys. 
Chem. B 109, 10542 (2005). 

[54] A. V. Ponomarev, J. Madronero, A. R. Kolovsky, and A. Buch- 
leitner, Atomic Current across an Optical Lattice, Phys. Rev. 
Lett. 96, 050404 (2006), |arXw:cond-mat/0509602| 

[55] L. Fedichkin, D. Solenov, and C. Tamon, Mixing and Decoher- 
ence in Continuous-Time Quantum Walks on Cycles, Quantum 
Inf. Comput. 6, 263 (2006), |arXiv:quant-ph/0509163| 

[56] P. E. Parris and P. Phillips, The role of dephasing in some re- 
cent theories of quantum localization, J. Chem. Phys. 88, 3561 
(1988). 



